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ABSTRACT 


,• 


Relationships  relating  power,  magnetic  field,  current  density,  and 
ampere-tunis  in  terms  of  certain  dimensionless  factors  are  suoimarized 
for  many  types  of  coil  geometries  and  current  distributions.  A  num¬ 
ber  of  plots  of  these  factors  are  presented. 

The  field  homogeneity  in  magnet  structures  Is  presented  in  lemos 
of  a  series  expansion  a’x>ut  the  origin  utilizing  Legendre  Polynomials. 

A  number  of  tables  to  facilitate  design  of  homogeneous  fields  are 
presented.  A  method  of  achieving  homogeneity  In  long  solenoid  struc¬ 
tures  by  the  use  of  determinants  is  discussed.  Expressions  for  the 
axial  field  from  uniform  and  radially  varying  .mrrest  density  coils  are 
given. 


« 
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PART  I 

RELATIONSIHPS  BETWEEN  MAGNETIC  FIELD,  POWER,  AMPERE -TURNS 

AND  CURRENT  DENSITY 


It  has  been  customaiy  in  the  literature  on  magnet  design  to  write  the  relatlonsulps 
between  fields,  power,  ampere-tums  and  current  density  in  tenna  of  a  number  o£  dimen¬ 
sionless  factors.  In  this  section  we  will  treat  in  summaiy  form  twelve  cases  of  coll 
construction  (Cases  I  -  XU).  This  work  is  In  part  a  summary,  and  enlargement  of  parts 
of  the  following  three  papers:  F.  Bitter,  R.S.  L  Vol.  7,  1936  Part  II;  W,  F.  Gauster, 
AIEE,  Fall  General  Meeting,  Chicago,  October  1959;  F.  Gaume,  Journal  de  Recherche 
du  Centre  Nationale  de  Recherche  Scientifique,  No.  43,  Jtme  1958.  (English  translation, 
Lincoln  Laboratory,  M81-12,  part  I  and  H,  Dr.  H.  H.  Kolm.) 


A.  The  most  commonly  used  dimenfaicnless  factor  is  the  G  factor  first  suggested 
by  Fabry,  which  connects  a  given  magnetic  field  with  the  power  required  to  produce  this 
field  and  is  a  function  of  the  type  of  current  distribution  and  the  geometry  of  the  coU- 


(1) 


with  W  =  watts 

p  -  resistivity  in  ohm  -  cm 

Gl  =  inner  radius  In  cm 

X  -  fractional  volume  of  conductor 

The  space  factor  A  is  assumed  coa.‘»tant  throughout  the  volume.  If  X  is  a  func¬ 
tion  of  X  and  y,  ^'ore  complicated  formulas  are  needed. 


D.  A  seco  V  (iimenslonlesa  factor  J  has  been  suggested  by  W.  F.  Gauster  and 
connects  the  curroni.  density  at  the  innermost  winding  of  a  mag^iet  to  the  power  of  the 
magnet.  J  is  •  function  of  the  current  distribution  and  the  geometry  of  the  coil. 


C.  i  v.  ....J  ..itlcular  type  of  coll  a  maximum  G  factor  can  be  found  which  will 
give  the  most  magnetic  field  for  the  least  power  used.  However,  !a  many  types  of  coll 
construction  this  maximizing  of  the  G  factor  may  lead  to  e-xcesslvely  high  current  densi¬ 
ties  at  the  innermost  winding  which  may  make  cooling  the  coll  difficult.  To  illustrate 
this  point  a  ratio,  X  •  forir.ed  which  compares  the  maxincum  current  density  In 

any  coll  to  the  current  density  In  a  uniform  current  density  coll,  where  both  have  the 
same  inner  radius  and  both  develop  the  same  field.  Quantities  for  fha  uniform  case  are 
designated  as  primed: 


,■  •♦r 


•t  amj^  •»'!  w  • 


r-'i  1  - ’-^^ 


H-»(®)“  H'-0{Sf 

if  H=  H' 


1  Q 

from  (2)  « 

J  {  wV/^ 
*  J'  \  w7 

and  using  (3) 

ji  "  r  G 


This  ratio  indicates  the  ditflculty  of  cooling  end  supporting  any  coll  compared  with 
a  coll  of  uniform  current  density  at  the  same  field  and  bore. 


D.  'fhe  ratio  in  equation  -(3)  Is  also  useful  and  represents  the  power  saving  in 
using  any  current  distribution  compared  with  uniform  current  density  when  generating 
the  same  field  in  the  same  bore.  We  give  It  in  the  more  useful  form. 

W'  \G  I  ^  (5) 

These  dimenslonlesn  factors  and  ratios  are  presented  in  Table  I  for  nine  types  of 
magnet  construction  and  following  Table  I  for  three  additional  more  complicated  magnet 
constructions. 

One  of  the  most  Interesting  new  coll  constructions  presented  is  that  of  case  V. 

This  coll  was  first  suggested  by  F.  Gaume  and  achieves  a  relatively  high  G  factor  with¬ 
out  excessive  current  densities  and  can  be  quite  simply  constructed  of  plates  in  the 
usual  Bitter  magnet  style,  with  the  introduction  of  axial  current  variations  by  making  the 
plates  thicker  towards  the  ends  of  the  magnet. 

A  number  of  plots  of  the  geometry  factor,  G,  against  CT  and  0  aie  given  In 
Figures  I  through  6.  Cockcroft's  curves  for  case  VlII  are  reprinted  In  Figure  1,  Bitter'a 
curves  for  case  VI  in  Figure  2,  new  curves  for  cases  vni  and  VI  in  the  small  CT  and 
0  region  are  given  in  Figures  3  and  4.  New  curves  for  case  V  (Gaume's  axial  current 
variation)  and  for  case  X  (the  ellipsoidal  shell)  are  given  In  Figures  5  and  6.  The 
variables  CT  and  j3  are  defined: 


I 


j  . 

( '• 

I- , 
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CL  =»  ^5^1 

%t 


E.  It  Is  sometimes  instructive  to  write  the  magnetic  field  in  terms  of  ampere- 
turns  rather  than  power.  This  rzises  from  the  practical  standpoint  that  It  is  difficult 
to  always  predict  the  resistance  of  a  magnet  exactly,  especially  as  a  function  of  power, 
and  one  may  wish  to  know  how  much  field  will  be  produced  for  each  ampere  available  In 
the  power  supply.  This  relatlorshlp  between  ampere-turns  and  field  can  be  related  to 
the  J  factor  and  the  G  factor  an!  miscellaneous  Ct  and  /3  terms  and  numerical  con¬ 
stants.  The  relationships  are  given  here  for  three  cases:  Case  VI,  Case  VII  and  Case 
Vin.  All  the  others  are  easily  derivable  by  Integrating  the  current  density  over  the  coil 
to  obtain  the  total  current  in  terms  of  J  and  W.  Total  current  I  and  the  factor  J  can  now 
be  substituted  for  W  in  equation  (1),  and  the  ampere-turns  related  to  the  field. 


Case  VI : 


H6 


Case  VII: 


Case  Vm: 
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?.  Useful  geometry  factorn  can  be  written  for  very  long  solenoids  (long  enough  so 
that  end  effects  can  be  neglected). 

For  a  constant  current  dennlty  coll  the  ampere-turns  per  unit  length  can  be 
written: 


Ni  /(w/2bn  a{£c-i)V/® 

zh^K-zH^J  i-5Tr7 

where  W/2b  is  the  power  per  unit  lenjjth  (watts/centimeter) 

For  the  disk  construction,  with  a  radial  current  distribution  we  write 


Rewriting  (9)  and  (10)  in  terms  of  the  field  In  gauss  and  power  per  length  as  watts 
per  meter,  we  obtain  .  .  ^ 


where 


11/  ^  /  W  (watts/rncter)  X 

G5  (uniform)  ^ 


G5  (radial)  =  ^J£m.  <3*)^^^ 


These  factors  are  plotted  in  Figure  7  and  it  is  interesting  to  note  that  below  an  CT 
of  2.  5  there  is  negligible  power  saving  with  the  disks,  but  at  greater  Ct  ,  the  saving  is 
Increasingly  pronounced;  at  an  CT  of  10,  45%  more  power  would  be  required  by  the  uni¬ 
form  wound  solenoid  for  the  same  magnetic  field. 

G.  It  Is  useful  to  use  the  G  factor  notation  when  dealL_g  with  split  colls  (colls 
separated  by  a  gap).  If  we  take  a  pair  of  coils  with  parameters  as  defined  in  Figure  8, 
we  can  write  the  field  In  the  gap  by  simple  addition  and  subtraction  of  coils. 


Figure  8 


(11) 


G: ^j)  f  G2  =Ga  (cTf 


This  relationship  (11)  holds  lor  any  current  dlatT-lbutlon  having  no  axial  variation. 
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CASE  X  {FIGURE  6) 

Coil  Shape;  Ellipsoidal  shell,  ^  =  constant 
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Geometry  Factor: 
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Gjnax  -204  at  cr  =  6,  /3  =  2 

Current  Density  Factor  J: 
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~  '  4  7r/i  in  cr 
at  cr-^,  >3=^2,  Gmctx  ^  0.204 
and  Jtjiax^OJ49 

Current  Density  Hallo  for  same  P'leid  referred  to  case  VIII.-  =  1.31 
I’owci  itatlo  P  for  same  1  leld  referred  to  case  Vlll  =  0.  763 
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CASE  XT 

Coll  Shape:  Cylindrical  Bore,  tapered  ends,  non-colccldlng  apexes. 


slope  “k 


b  =  ky  -f-o 
k  =  slope 
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Current  Distribution: 


Geometry  Factor: 


i  -  io/y 


1 _ _  r  k 

5  l<(cr-0+ 


fl  ((k"+l')a'^+2kya'  +  y2)'/2  +  a’  (k^+l)  +  y 
((k^+1)  +  2ky +y2)‘/2  4  (k^fl)'/^  f  + 

^Jui or  4  y  4  k _ ^“1 

\  (Ck^fl)  a:^42Jcycr4y2)‘''^4y 


•  9 


9  9 


9  9 


G  is  a  maximum  when  k  =  1  and  CC  =4.5  which  reduces  case  XI  to  case  VH  Current 


Density  Factor  J: 
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CASE  xn 

Coll  Shape;  Cylindrical  Bore,  tapered  ends,  non-colncldlng  apexes. 
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Current  Distribution: 


Geometry  Factor; 
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G  Is  a  maximum  when  k  =  1.0  and  CC  =4.5  which  reduces  case  XII  to  case  VII  Current 
Density  Factor  J: 
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PART  n 

HOMOGENEOUS  MAGNETIC  FIELDS 

A. 

A  magnetic  field  cr  n  be  written  as  a  power  series  involving  I,egendre  polynomials 
expanded  about  the  origin  of  the  field.  If  this  origin  is  on  the  axis  and  on  the  plane  of 
symmetry  the  expansion  will  have  no  odd  terms  and  the  axial  component  Hz  and  the 
radial  component  Hyo  can  be  written: 


H:=  -I 

n-=  i  L  J 


2n-2 


2n-2 


H/3  ^  Z  (2n_1)i  U|0)]  ^  P2r»-2  <2) 

U  =  COS  ^ 

In  these  eqiiationa  the  coefficients  are  defined  as  follows: 


and  p2n-2(*^)^^®  Legendre  polynomials  tabulated  in  Table  2. 

The  importance  of  this  approach,  due  to  W.  Garrett^,  la  that  the  field  at  any  point 
on  the  axis  of  a  system  of  coaxial  colls  can  be  given  explicitly  in  simple  algebraic  form. 
From  this.,  the  derivatives  of  equations  (1)  and  (2)  can  bo  obtained,  and  so  a  complete 
expression  for  the  field  in  the  vicinity  of  the  origin  can  bo  obtained.  If  sufficient  deriv¬ 
ative  terms  are  used  the  field  can  bo  found  an3rwhore  within  a  sphere  whooe  center  lies 

■  t 

at  the  origin  and  which  does  not  extend  far  enough  to  Include  any  comers  of  the  coil.  How¬ 
ever,  out  to  within  a  few  percent  of  the  Imjer  radius  of  the  coil,  the  field  can  be  fouud 
quite  accurately  with  only  a  few  derivative  terms.  Compensation  can  be  designed  for  each 
component  of  the  field  separately.  In  general,  if  wo  have  a  coll  system  of  the  assumed 
symmetry  having  N  variable  parameters  such  as  size,  or  ampere-turns.  It  Is  possible  to 
specify  N  coefficients  in  the  expansions  1  and  2. 

For  example,  to  examine  the  standaid  Helmholtz  pair  of  colls,  where  the  colls  are 
so  spaced  that  there  is  no  second  derivative  of  the  field  at  the  origin,  we  can  think  of 
the  coil  system  In  tlie  following  way: 


1. 


M.  W.  Garrett,  J.  Appl.  Phys.  22,  9,  Sept.  1951 
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Co)  Cb)  tc) 

C*.  b  “G 

Figure  1 

As  shown  in  Fig.  1,  v;e  consider  a  single  coll  having  length  equal  to  the  coils  and  the 
gap  separation,  from  which  we  svibtract  a  second  coil  of  the  length  of  the  gap.  The  gap 
coil  is  thought  of  as  having  a  current  in  the  reverse  direction  to  the  original  long  coil. 
We  now  merely  need  select  the  length  of  the  gap  coil  so  that  it  will  have  a  second  deriv¬ 
ative  equal  and  opposite  to  that  of  the  original  long  coll  and  the  two  will  thus  cancel  at 
the  origin. 

There  are  many  other  combiaatlona  which  would  also  cancel  the  second  derivative 
at  the  origin,  for  inaUuice,  the  hypothetical  gap  coil  with  its  reverse  current  could  have 
been  somewhat  longer  and  had  less  reverse  current  in  it,  or  shorter  and  had  more  re¬ 
verse  current  In  it.  Both  of  these  solutions,  of  course,  lead  to  a  continuous  solenoid 
with  no  gap  but  with  an  altered  current  density  in  the  center  section-  These  various 
solutions  will  have  both  an  effect  on  the  amount  of  field  in  the  center  of  the  magnet  and 
on  the  size  of  the  higher  derivatives.  The  greater  the  distance  over  which  one  compen¬ 
sates,  the  more  parallel  will  be  the  flux  lines  from  the  compensating  coil,  hence  the 
flatter  will  be  the  flux  distribution  at  the  origin,  and  hence  the  smaller  will  be  the  sum  of 

the  higher  derivatives.  This  compensation  may  not  be  compatible  with  obtaining  the 

2 

highest  central  field  for  the  least  power,  however. 

B. 

Equations  (1)  and  (2)  can  be  rewritten  in  the  following  form.-: 

H;.  (r,o)  =.  Hz  (0,0)  [l  r (u)  +  £4  (li)  P4  (u)  +  -]  <^> 

(n,o)  *  Hz  (0,0)  [^Ot£z  (^^^2  (“)+£4  (;^)  P4  (“)+  " '  j 


2-  R.  S.  Ingarden  and  J.  Mlchalczyk,  Bulletin  De  L’Academie  Polonaise  Des  Sciences 
(Scrle  des  scl.  math,  astr. ,  and  phys)  VIU,  5,  1960. 
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The  expansion  has  been  normalized  to  the  Inside  radliis  of  the  coils  A  tabulation  of 

the  coefficients  for  three  useful  geometries  is  given  In  Table  1.  The  geometries 
treated  are  (1)  a  cylindrical  thin  current  sheet,  (11)  a  finite  solenoid  with  a  unifonn. 
current  density,  and  (iii>  a  finite  solenoid  with  a  radial  current  distribution  (3^  %>. 

The  value  of  the  first  eight  Legendre  polynomials  and  their  derivatives  are  given  for  con¬ 
venience  in  Tal'le  2,  From  equations  (3)  and  (4)  the  magnetic  field  can  be  found  at  any 
point  near  the  origin  of  a  magnet,  at  a  distance  V  from  the  origin  and  at  an  angle  ^  . 
The  derivative  terms  in  equation  (3)  have  their  largest  value  for  =  0  (P^  (1)  =1  for  all 
n  when  -  0  and  cos  ^  - 1).  The  maximum  value  of  the  terms  therefore  can  be  ob¬ 
tained  from  the  simpler  form  of  (3)  where  ^  -0  and  Z  =  P - 

Hz  (Z,o)  =  Hz  (0|0)  Ij  ^^2.  (^,)  ^  ^4-  +  "'j  (5) 

c. 

The  most  oominonly  used  method  of  achieving  some  degree  of  homogeneity  in  colls 
is  to  use  the  spaced  heircholtz  pair.  Garrett  has  compiled  a  number  of  tables  dealing 
with  this  method  of  corr.psnsation  and  since  these  tables  have  not  been  published  we  are 
making  them  part  of  this  report.  These  tables  are  for  constant  current  density  only. 

The  first  table.  Table  3  is  used  to  determine  the  proper  spacing  between  finite  colls  to 
achieve  cancellation  of  ths  second  derivative  at  the  origin.  The  second  table,  Table  4 
Is  used  to  determine  honr  far  from  the  origin  one  can  move  before  the  field  deviates  by 
more  than  0. 1  percent  of  the  field  at  the  origin.  The  third  table,  Table  5  is  particularly 
useful  and  can  show  how  accurately  the  colls  must  be  spaced  in  order  to  cancel  the 
second  derivative  to  the  desired  degree.  The  two  remaining  tal)les,  Tabie  C  and  Table  7 
give  the  resultant  fourth  ar.!  the  sixth  coefficients  x'emaining  after  the  proper  spacing 
of  the  coils.  The  tables  are  used  with  the  notation  of  equations  (3)  and  (4)  except  that 
Garrett  has  normalized  ths  exfxmslon  to  the  mean  radius  of  the  coll  rather  than  to  the 
inner  radius  as  in  Section  B  above  (  (-^)  rather  than  (.^  )  whore  <3^  =  ). 

To  illustrate  the  use  of  G.irreit's  tables,  we  work  out  the  following  examples: 


•  • 


•  • 


Assume  we  have  tv.o  coils,  each  with  an  CT  -  3  and  n  1  and  wo  wl.sh  to  sepa¬ 
rate  thorn  to  cancel  the  swond  order  derivative  at  the  origin. 


X  - 
A- 


4/3  ^ 

a*+i 

2(0-1) 

a'  + 1 


1.0 

-  1.0 


•  • 


•  • 


From  Table  3, 


K  -&^2ro 


AX 

Clj  +  Q^ 


AX, 


•  • 
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therefore 


Z.33040  «  2y3  4-  ^ 

‘  ai 


and  the  spacing  should  be 

=.  .;.3040 

ai 

Using  Table  4  we  see  thatV=^  =181  or  ~  •  =.724.  This  means  ihat 

out  to  a  distance  along  the  axis  of  72.4%  of  the  Inner  radius,  the  field  Is  within  0.1%  of 
the  -field  at  the  origin.  Using  Table  5,  we  can  examine  the  effect  of  Improper  spacing. 


where  Q  Is  the  value  In  table  5  and  d/a^  Is  the  displacement  error  In  locating  one  coil. 
For  the  above  example,  -we  can  find  the  error  that  will  make  the  uncancelled  2nd  deriva- 
tlve  error  equal  the  fourth  at  (  ~  )  =  0. 1 

Q  =  4.21  (at  X“1,A  =  l) 

<$^4=  -922  from  table  6 
4.2t  X  10"^  ==  -922  X  10“^ 

^  =  2.18  X  10"^ 

ao 

D. 

When  maximum  homogeneity  and  not  maximum  efficiency  la  of  the  greatest  impor- 
Lance,  the  higher  derivatives  must  be  handled.  Since  the  number  of  derivatives  that  can 
be  cancelled  depends  upon  the  nuniber  of  variables,  we  must  now  pick  solutions  which 
have  more  variables  at  our  disposal.  For  Instance,  If  wc  wish  to  c.ancel  both  the  second 
and  the  fourth  derivative  of  the  field  V'C  can  GujKirlmpo.so  two  coll.s  on  each  other  each 
with  two  geometry  variables  or  each  with  one  geomotiy  variable  and  one  current  variable. 
To  proceed,  then, we  would  pick  the  Icivjth  of  one  coil  and  the  current  of  one  coll  and  the 
variables  then  v/ould  bo  the  length  of  the  second  coll  and  the  ratio  of  the  two  current 
den.sltles.  Garrett  has  calculated  a  set  of  so  called  sixth-order  solenoids  (no  second  or 
fourth  derivatives)  by  inaKlng  Uic  following  choice  of  variables:  He  constructs  an  over¬ 
wound  end  solenoid  by  winding  N  layers  of  turns  over  a  certain  central  region  of  the  sol¬ 
enoid  and  2N  layers  over  a  certain  section  of  both  t-ntis  (l.e.  a  notch  on  the  center 
section  of  the  o.  d.  of  the  coll).  Ho  now  .solvo.s  uniquely  for  Uio  length  of  the  overall 
coll  ami  the  lonf^^h  of  the  notch.  The  results  of  his  calculation  are  shown  In  Table  8  for 
several  t.'ilcknes.sos  of  coll. 

There  are  ottier  cholce.s  of  two  variables  that  ono  cun  laako  to  con.struct  sixth-orxitm 
suleiiold.s  and  Table  0  Illustrates  se^'oral  inethods.  T.iblo  9  was  constructcii  for  the 
.specific  example  of  a  ra'her  b'ug  solenoid.  Several  gcn-'rul  ren^arks  can  be  dr.awn  from 
the  table.  I’robaldy  the  most  liTij-oifant  of  which  is  that  since  the  solenoid  is  very  long 


•  • 


the  job  of  ccmpenaatlng  It  for  high  homogeneities  Is  much  simplified,  the  higher  derivative 
having  been  quite  small  to  begin  with.  We  notice,  however,  that  In  all  but  one  case  In 
the  table,  cancelling  the  second  and  fouilh  derivatives  of  the  field,  have  Increased  the 
sixth  derivative.  The  single  exception  l,o  the  current  sheet  of  the  o.  d.  The  current  sheet 
on  the  o.d.,  however,  requires  more  power  than  any  of  the  ether  methods  of  solution. 
Perhaps  the  most  attractive  method  of  solution  Is  that  represented  by  the  last  column, 
that  of  current  sheets  on  either  end  of  the  l.d.  of  the  coll.  This  method  takes  little 
power  and  leaves  the  center  section  of  the  magnet  open  for  further  compensation  If  desired. 

Cancellation  of  higher  derivatives  than  the  fourth  la  of  course  possible,  but  la  often 
of  decreasing  practical  Importance  because  of  the  unrealistic  accuracy  necessary  In  the 
cancellation  of  lower  derivatives,  If  the  leftover  uncompensated  part  of  say  the  second.  Is 
to  be  smaller  than  the  sixth  derivative  or  eighth  derivative,  etc.  It  probably  makes  sense 
to  design  colls  to  cancel  the  second  and  the  fourth  derivative  and  use  a  small  auxiliary 
set  of  colls  driven  from  the  separate  power  source  to  tuno  out  any  random  errors  which 
creep  In.  Some  of  these  random  errors  are  of  course  going  to  be  nonsymmetrlcal  and 
will  give  rise  to  odd  derivative  terms.  It  la  quite  possible  to  construct  a  set  of  auxiliary 
compensating  coils  which  will  operate  Independently  of  each  other.  One  coll  can  be 
designed  to  have  a  second  derivative  and  no  fourth  and  a  second  coll  to  have  a  fourth 
derivative  and  no  second,  etc.  Each  one  can  then  be  tuned  Independently  of  the  other. 

E. 

We  have  constructed  a  table  of  €  coefficients  up  to  the  eighth  order  for  coils  of 
varying  O’  and  /i  with  both  constant  current  density  and  with  current  density  equal 
i  '‘Vr  .  The  tables  also  Include  the  G  factor  and  the  multiplication  of  the  G  factor 
and  the  &  coefficients.  The  tables  are  to  be  used  with  the  expansion  of  equation  (3) 
and  (4)  where  z  Is  referred  to  the  inner  radius  Cl j ;  <5  Is  tabulated  as  D  for  uniform 
currents  a’ul  as  A  for  the  radial  current  distribution.  Table  10.^  gives  the  coefficients 
for  the  uniform  current  demslty  ca.se  for  several  OT  3  and  for  up  to  10.  Table  lOB 
give.s  the  coefficients  over  Uie  same  range  of  CC  and  /3  for  the  radial  current  distribu¬ 
tion  case.  Tables  llA,  B,  C,  D  show  more  detail  for  the  constant  current  and  the 
'•adlal  current  cases  in  the  small  range. 

The  use  of  the  tables  Is  probably  beat  Illustrated  by  the  following  example:  I,et 
u.H  .assume  wo  have  a  coll  we  wish  to  separate  for  pui]x>r;es  of  access  .and  we  want  to 
know  the  resultant  f^old  In  the  g.ap  and  the  resultant  homogeneity  (see  Fig.  2) 
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X  =  space  factor 
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Figure  Z 
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If  v,e  now  subtract  from  the  field  expansions  for  a 
i  we  can  write: 

coll  of  length  4b  +  ,  a  coll  of  length 

9  9 

^gop  ~  (^ot” 

^^0g)+(H2t  ”*^2g) 

(^)  ■»'  (^4t  '  +-  "* 

where  the  primed  qmiirtules  refer  to  the  error  coefficients  for  the  gap.  TaVlng  a  spe¬ 
cific  geometry,  let  the  current  he  corstant  {&=  D). 
let 

/3t  =  4 
O'  ='  3 

/?g  —  0.2 

Whzit  is  the  field  in  the  gap  and  the  second  error  coefficient,  and  what  would  be  the 
second  error  coefficient  If  there  were  no  access  gap? 
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•  • 


•  • 
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.-1 


The  Di  listed  for  OT  =  3,  /3  -  4  In  Table  lOA  is  -  .1415  x  10  *,  therefore  the  second 


c*  -2 

error  coefficient  of  the  whole  coll  would  be  02=  -1.415  x  10  where 

H-'Ho  (•  +  "') 

To  solve  for  the  split  coil,  we  use  (6)  and  the  G| .  Dn  columns  in  Table  lOA 

-8- 0.2235  X  10"^+  0.5004  x  10“^  '"J 

“(^n(*  05)[  <387 4-0.8965 

this  gives  a  central  field  of 

Hgap  (0|0)  -  (o.f43) 

and  If  the  second  order  coefficleEt  Is  written  in  the  following  form: 

H  =  Uo  (1  +  £z  (|)%  —) 

the  coefficient  for  the  above  example  would  be 

=  5.46  X  10“^ 

A  387 

The  loss  ill  field  for  this  case  has  been  9.  5%  and  the  second  order  error  coefficient  has 
been  increased  by  a  factor  of  4.55,  If  we  had  been  trying  to  homogenize  the  field 
separating  the  coils,  we  note  that  we  would  have  separated  them  by  too  much.  To  "le 
these  tables  for  synthesis  (1.  e-  cancellation  of  derivatives)  it  would  be  necessary  to  use 
trial  and  error.  In  the  manner  of  the  above  example. 

It  Is  interesting  to  note  the  oscillatory  behavior  of  the  <5  coefficients,  particularly 
lu  the  small  y3  region.  The  second  derivative  Is  always  negative,  the  fourth  derivative 
goes  through  zero  once,  the  sixth  derivative  goes  through  zero  twice  and  the  eighth 
derivative  three  times.  This  meatm  that  there  are  colls  which  can  be  built  with  no  fourth, 
no  sixth  and  no  eighth  derivatives.  This  may  bo  useful  In  the  design  of  Independent  com¬ 
pensating  colls. 


Sma  .i  ‘liianiTM  ■TI’A'i*.  .fawairL'  . 


F. 

Iff  one  wishes  to  compensate  a  very  long  solenoid  over  a  considerable  lerigth  of  the 
axis  such  as  a  multicoil  structure  used  for  plasma  experiments,  expansion  of  the  deriva¬ 
tives  aroimd  the  origin  has  limited  usefulness.  One  might  wish  to  use  the  following  pro¬ 
cedure  instead;  namely  specifying  that  the  fleid  at  any  number  of  intervals  along  the  axis 
be  constant  and  then  solving  for  ti  e  current  distribution  in  the  colls  that  would  make  the 
field  be  equal  at  the  specified  points.  In  the  case  of  a  multlcoll  structure  one  could  ask 
that  the  total  field  in  the  center  of  each  of  the  coils  be  equal  to  a  constant  and  then  solve 
for  the  current  density  required  in  each  of  the  colls  to  make  this  true.  The  same  proce¬ 
dure  could  be  used  with  a  single  loug  solenoid  where  the  solenoid  could  be  broken  up  into 
an  arbitrary  number  of  pieces. 


The  most  useful  method  of  solution  of  this  type  of  problem  Is  basically  the  following: 
taking  any  single  ceil  we  first  find  the  field  at  the  center  of  that  coil  and  at  distances 
away  from  the  coll  represented  by  the  center  lines  of  each  of  the  other  coils,  for  a  given 
current.  Then  by  properly  adding  up  the  coils  In  matrix  form  and  demanding  that  the 
field  matrix  times  the  current  matrix  be  a  constant  we  can  invert  the  current  matrix  and 
solve  for  the  individual  currents  required  In  each  of  the  colls  to  achieve  the  desired  re¬ 
sults.  The  method  of  solution  and  a  useful  notation  are  developed  in  the  following  example. 
Mr.  R.  Bradshaw  of  Conesco,  Inc. ,  Arlington  ,  brought  this  method  of  solution  to  our 
attention.  Consider  the  following  system  of  coils  arranged  so  that  they  have  a  common 
axis.  They  may  in  general  be  of  various  sizes  but  for  Illustration  consider  them  to  be 
all  identical  with  the  same  separation.  For  ease  in  calculation  choose  an  even  number 
of  colls  so  that  we  have  a  plane  of  symmetr'/  passing  between  colls. 


in  each  coil 


Our  object  is  to  find  what  ratio  of  currents  is  needed  In  the  various  coils  so  that  a  uni¬ 
form  magnetic  fie*  ’  ja  obtained  throughout  the  interior.  If  the  coils  are  thin  and  close 
together,  to  i  ;v  approximation  uniformity  is  obtained  by  requiring  that  the  field  at 
the  center  of  e  coll  be  the  same.  The  field  falls  off  rapidly  on  each  side  of  the 
center  of  a  s  '  ■  and  we  may  construct  a  matrix  the  elements  of  which  represent 
the  contrlbutiona  o.  I  !o  coils  at  their  respective  centers.  Lot  Fij  be  the  field  at  the 
center  of  coll  t  duo  to  the  field  of  coll  j  .  For  exan»ple,  for  4  coils,  2  colls  on  each 


•  • 


side  of  the  plane  of  symmetry,  we  have 

\ 

^14  ^12  •  fl3  Fi4 

Fei  F22IF23 

F32  j  F33  F34 

I 

F41  ?42  I  ^43  ^4 

I 

These  matrix  elements  may  be  expressed  in  %  of  an  Infinite  solcnoldal  sheet  of  the  same 
mean  radius,  or  whatever  one  desires,  since  we  need  only  ratios. 

FiJ  Ij  =  Tx  where'll  Is  the  field  at  the  center  of  the  i^^coil  or 

I 

Ffi  Ft2  1  Fi3  Fi4  It  y, 

I 

F2I  F22  I  F23  F24  Iz 

F32  J  F35  F34  I3 

Hta  I  F43  Kf4  I4 

If  one  now  specifies  *ae  ^1*3  the  Fjj  matrix  can  be  Inverted  and  the  Ij  matrix 

found. 

If  the  current  in  the  coils  are  arranged  so  that  It“l4  arid  I^, "I3  then  It  is  convenient 
to  partition  the  matrix  as  follows; 
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f  _  hi  F'Z  r  _  F,5  F,4  F3,  P32  f  Fjj  Fm 

~  F.,F22  ■  F23F24  ’  F4,F«  •  fe-  F43F44 


JFjjFm 


•  • 
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•  • 
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i,  -  ^3  , 


^  r7i 


•  • 


fll  fl2 

- 

and  then  '2  - 

fez 

5^2! 

where  k  »j 

0  1  Ii 


Til 

fl2 

tl 

^21 

(22 

kis. 

•  • 


it  follows  that 


(fn  +  fi2k)i,  =  y; 


r- 


-.ii 


jfvwf ’ 


Since  we  desire  the  fields  at  the  centers  of  each  coil  to  be  a  constant 
“  3^2.  ~  constant  =  C  whore  C  C  |  |  j 

i.e.,  (fn  +  ^iKi^)ii  ~C 

It  Is  clear  that  i|  and  this  gives  us  the  desired  ratio  of  currents. 


As  an  example  consider  the  field  on  the  axis  due  to  a  system  of  four  current  loops 
of  radius  y  =  1  Cm 


11-  I4 

12- I3 


Figure  4 


Separation  between  centers  of  the  coils  is  1  cm.  We  need  calculate  the  field  of  one 
loop  at  the  center  and  1,  2,  and  3  cm.  on  the  axis  away  from  the  center.  For  a  single 
loop 


then 
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27ri 

10  (x2Hy2)3/2 


0,620  0.222 

0.222  0.628 


fl2k  = 


0.056  0.019 
0.222  0.056 
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X 
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0 
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1 

.222 

2 

.0.562 

3 

.0199 
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0.222  0.056' 


0.019  0.056 
0.056  0.222 


0.628  0.22Z 

0.019  0.056 

0.647 

0.278 

0.222  0.628 

0.056  0.222 

0.270 

0.850 

We  now  need  the  lnver.se  of  (  ^'ll'1'^12^) 


By  definition,  tha  Inverse  of  a  matrix  A  with  matrbx  elements  ajj  Is  obtained  If  for 

each  element  (Hij)  ^  j  Aji  is  the  minor  of  the  determinant  A. 

det  A 
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a  A  ■■■I'*"  then  A"’- 

1  021  «22i 


For  2  X  2  matrices  this  Is  very  simple: 


(oit)  \atz)  * 

(<^21)  ^C^zz)  ^ 


A-'- 


det 


\ _ /  ^zz  ■* 

tA  \~a2i  an/ 


det  (fii  +  ftzk)  -  (.850)(.647)  -(27S)2  =*  .473 


.850  -.278 

-.278  .647 

1.797 -.587 
-587  1.367 


! 

1 

If 

1.797 

-587 

^ 

1.210 

•  • 

; 

l2 

"  -.587 

1.367 

1  c  c 

.780 

A  Ii  ^  C  (1.210)  l2^C  (0.730) 


The  current  In  coil  1  must  then  be  1. 55  times  greater  than  the  current  In  coll  2. 

Ingarden  and  Mlchalczyk^^^  have  also  considered  the  problem  of  homogeneities  over  a 
certain  length  of  the  axis  rather  than  strictly  at  the  origin.  They  have  used  the  require¬ 
ment  that  the  mean  square  deviation  over  a  certain  inten'al  be  minimized.  Their  notation 
is  useful,  but  solution  of  the  equations,  which  are  transcendental  In  character,  presents 
a  number  of  problems, 

G. 

It  Is  often  of  Interest  to  know  the  field  along  the  ajtls  of  a  solenoid  both  Inside  the 
solenoid  and  outside  the  solenoid.  'Fhe  following  two  equations  can  be  used  to  find  the 
field  along  the  axis  of  a  uniform  current  solenoid  and  a  radially  distributed  current  sole¬ 
noid  (Bitter  discs)  respectively. 

(7) 

kZM/2 


(i)  Uniform  current  % 


TTiXaiffl 
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(ii)  Radial  current  ^  ^ 


ai 


H  -  ^  (k+a)  I ,  t (k+/))2)'«* 


(S) 


Tabiilated  values  of  equation  (7)  for  a  wide  variety  of  coll  shapes  have  been  published  by 
D.  E.  Mapother  and  James  Snyder  In  a  University  of  Illinois  circular  #66  entitled,  ’’The 
Axial  Variation  of  the  Magnetic  Field  in  Solenoids  of  Finite  ITilckness”.  There  la  no 
available  tabulation  of  equation  (8) 

Equations  (7)  and  (8)  can  be  written  In  a  much  simpler  form  utilirlng.  the  following 
angle  notation: 
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I)  uniform  current  density 

TriiXai 

H  =  — _ 


+ 

4  a  / 


(9) 


ii)  radial  current  density  p 


1 1  _  vri'tXat 


tan  ^  tan  ^ 


tan  - 


^ii 


fan 


(10> 
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General  expressions  for  the  off  axis  fields  arising  from  a  loop  of  current  can  be 
written  in  terms  of  elliptic  integrals  of  the  first  and  second  kind.^^^  They  can  be  written 
in  the  following  way,  where  the  parameters  are  defined  In  Figure  8. 
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(11) 


(12) 
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W.  R.  Smythe,  STATIC  AflD  DYNAMIC  ELECTRICITY,  p.  266,  McGraw-Hill,  New 
York,  1950 
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Elliptic  integral  of 
the  first  kind 
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where  E(k)  = 
elliptic  integral 
of  the  second 
kind 
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Using  equations  (11)  and  (12)  the  field  at  any  point  produced  by  a  finite  coll  can  be  found 
by  dividing  up  the  coll  Into  loops  and  summing  all  the  contributions.  A  number  of  labora¬ 
tories,  including  our  own  have  written  simple  computer  programs  utilizing  equations  (11) 
and  (12)  to  find  the  fields  in  and  outside  the  conductor  volumn. 


On-axis  and  off-axis  fields  from  any  solenoid  with  any  radial  current  distribution 
can  be  calculated  by  means  of  a  lengthy  set  of  tables  published  by  the  Oak  Ridge  National 
Laboratory  Report  #ORNL  2828,  entitled,  "Tables  for  a  Semi-Infinite  Circular  Current 
Sheet,  "  by  N.  B,  Alexander  and  A.  C.  Downing.  These  tables  are  for  the  field  around 
the  end  of  a  seml-inflnlte  current  sheet  and  the  coll  in  question  is  built  up  out  of  these 
current  sheets  in  the  method  shown  in  the  introduction  of  the  table. 


Attention  Is  also  called  to  a  tabulation  by  E.  E.  Callaghan  and  S.  H.  Maslen,  NASA 
D465  "The  Magnetic  Field  of  a  Finite  Solenoid".  This  paper  plots  graphs  of  the  radial 
and  axial  components  in  and  aro'ond  current  sheets  of  various  lengths.  Wlille  not  as 
accurate  as  the  Oak  Ridge  tables,  it  is  very  useful  for  approximate  results. 


A  final  paper  of  particular  relevance  Is  that  of  G.  R.  North,  "Some  Parameters  of 
Lumped  .Solenoids",  (Oak  Ridge  ORNL-2975).  He  derives  an  expression  for  the  field 
ripple  resulting  from  the  separation  of  coils  in  a  long  msiltlcoll  structure  (end  effects 
neglected).  The  equation  Is  as  follows  (.see  Figure  7) 


3tn 


2b 


2bir 


•  cos 


Zttz 


2frr' 

'~s~ 


] 


l/2  - 

cr  e 


2-rra:a\ 


} 


(13) 


I 


TABLE  1 

TABLE  OF  ERROR  COEFFICIENTS 
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'VABLE  1  (continued) 
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